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Abstract. 

We present a study of the motion of compact jet components in quasar B3 1633+382. Through analyzing 14 epochs of VLBI 
observations of three components (Bl, B2, and B3) at 22 GHz, we find two different possibilities of component classification. 
Thus two corresponding kinematical models can be adopted to explain the evolutionary track of components. One is a linear 
motion, while another is a helical model. Future observations are needed to provide new kinematical constraints for the motion 
of these components in this source. 
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1. Introduction 

Associated with a redshift of z = 1.814 
dStrittmatter et al. 1974 . the radio s ource B3 1633+382 (4C 
38.41) was identified as a quasar JPauUnv-Toth et al. 1973 ) 
and has been frequently observed at different wavelengths 
across the electromagnetic spectrum during the past decades. 
A nearly flat spectrum has been deriv ed from 1.5 t o 90 
GHz contemporaneous measurements (ILandau et al. 19 86). 
Strong variability at both long and short radio wave- 



lengths has als o been detected dSpangler & Cotton 1981. 
lAlleretal. 19"92l: ISeielstad et al. 19851: iKuhretal. 198lb . 



This quasar has been classified a s an optically violent 
variable (OVV) (iMattox et al. 1993h because of its strong 
optical variability. Both infrared and X-ray fluxes have been 
obtained w ith the IRAS and the Einstein Obser vatory re- 
spectively dNeugebauer et al. 19861: iKuetal. 19801) . As one 
of the most powerful extragalactic objects detected by the 
Energetic Gamma Ray Experiment Teles cope (EGRET) on 
the Compton Gamma Ray Observatory ( Fichtel et al. 1994 : 



Thompson et al. 1995t iHartman et al. 19991) . its y-ray power 



is two orders of magnitud e higher than the typical value at 
any longer wavelength (Matto x et al. 1993b . This source is 
also included in the Fermi-LAT 3-month bright AGN list 
dAbdo et al. 20091) . 

Very Large Array (VLA) observations of B3 1633+382 
have shown a triple structure with an overall size of 12 
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along the north-south direction. This is a misaligned source, 
with its VLA jet at almost 90 degrees p osition angle to the 



Very Long Baselin e Array (VLBA) jet (IMurphv et al. 1993 



Britzen et al. 20071) . A jet in the western direction with a 
sha rp bend to the sou th at about 50 mas was detected 



by (IPoIatidis et al. 19951) with 1.7 GHz Very Long BaseUne 
Interferometer (VLBI) observations. Superluminal proper mo- 
tion jj. - 0.16 + 0.03 mas yr"' (4.4+0.8 c) was detected 
on th e basis of three e pochs of VLBI observations at 5 
GHz (Barthel et al. 1995). It has also been observed with the 



VLBA at dual-frequency (Fev & Chariot 1997h . The VLBI 
space observatory programme (VSOP) mission, using the 
HALCA satellite along with a global array of earth-based tele- 
scope s, has detected a faint core with su b -milliarcsec r esolu- 
tion CListeretal. 20011 IScott et al. 2003) . Ijorstad et al .1 (2001) 
presented results of their multi-frequency monitoring program 
with the VLBA during six epochs between 1994 and 1996. A 

prominent jet at a position angle (P. A.) 87°, extending out 

to 1 .5 mas with three components, has been detected at 22 GHz. 
The authors found that the proper motions of components B 1 
and B3 are 0.20 and 0.14 mas yr ' (5.5 c and 3.9 c), respec- 
tively. Component B2 was shown to be quasi-stationary at a 
distance of ~ 0.5 mas from the core. 

In our 22 GHz VLBA observation preformed on 2000 
March 1 (see Fig. 1), there are four jet components (in this 
wor k, we name them components Bl, B2, and B3 follow- 
ing Ijorstad et al. 200 ll) loacted in a western direction with 



0.51, 1.33 and 2.15 mas separation from the core respectively. 
Component B4 is the newest one at 0.15 mas from the core. 
This result is co nsistent with another VLBA observation at 
43 GHz in 1999 dLister & Smith 2000l) . However, after com- 
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bining their proper motion with Jorstad et al.l f2001') and ours 
(see Fig. 2(a)), we find it hard to distinguish the evolution- 
ary track of the two internal components owing to the lack of 
observations from 1997 to 2000. In order to study the kine- 
matic in quasar B3 1633-1-382, more observation epochs are 
needed. Thus we additionally collected and analyzed seven 
more VLBA epochs at the same 22 GHz frequency from the 
National Radio Astronomy Observatory (NRAO) archive data 
(see Table 1). Combined with these seven archive observed 
data, one observation of our own, and six other observations 
obtained from literature, we have 14 epochs of VLBA observa- 
tions at 22 GHz for this source. 

We study the kinematics of the compact jet in quasar B3 
1633H-382. The paper is organized as follows. In Sect. 2 the 
data reduction is described. The motions of components are 
detailed in Sect. 3, followed by our conclusions in Sect. 4. We 
adopt the spectral index convention /v oc y^" and a cosmology 
with Hq = 70 km .r' Mpc-\ Q.m = 0.3, = 0.7. 



I 
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Fig.l. Model-fitted image of the quasar B3 1633-1-382 at 
22 GHz, observed on 2000 March 1. The restoring beam is 

0. 737x0.319 mas at RA. = -1.6°. The contour levels are (-1, 

1, 2, 4, 8, 16, 32, 64, 128, 256) x 2.93 mJy beam-\ The peak 
brightness is 0.992 Jy beanT^ . 



2. Data reduction 

Our 22 GHz VLBA observation of B3 1633-H382 was per- 
formed on 2000 March 1. The signals were recorded in four 
intermediate frequency (IF) bands for a total bandwidth of 32 
MHz with 2 bit sampling. The recorded data were first cor- 
related at the VLBA correlator in Socorro (New Mexico) and 
then calibrated and fringe-fitted using the NRAO astronomi- 
cal image processing system (AIPS) software. Initial amplitude 
calibration was done using system temperature measurements 
and the NRAO-supplied gain curves. The data were corrected 
for 2 bit sampling errors and atmospheric opacity. The imag- 
ing and self-calibration w ere carried out in the DIFMAP pack- 
age (Shepherd, Pearson. dShepherd et al. 1994 . Of the seven 
epochs from the NRAO archived data, four were observed at 
8, 15, 22 and 43 GHz, one at both 22 and 43 GHz, and two at 
22 GHz only. All archive data were processed with the steps 



described above. However, we encountered a problem with the 
correction for atmospheric opacity for the data of 1998 January 
2. We used the simultaneously observed 15 GHz data instead 
to replace the 22 GHz data in our study, which might introduce 
a less than 0. 1 mas position shift. 
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Fig.l. Core separation as a function of time for the individ- 
ual jet component. The data combined are from the Jorstad et 
al. (2001) (open symbols), archive data, and our observation 
(solid symbols). Stars show the components from the VSOP 
observations by Lister et al. (2001). The first panel plots the 
proper motion with Jorstad et al. (2001) and our own obser- 
vation (purple pentagon). The second shows a simple linear 
proper motion to fit the evolution track of components of 14 
epochs. The last one displays a simple helical fit to the evolu- 
tion track of components depending on a different classification 
of the components. 

After obtaining images in DIFMAP, we performed the 
model fitting to the self-calibrated visibility data using circu- 
lar Gaussian components. For all model fits we assumed the 
brightest component to be stationary and used it as a refer- 
ence point. The other components were measured relative to it. 
Uniform weighting of the uv data was used to ensure the high- 
est resolution. Figure 1 shows the naturally weighted image of 
our observation with the model-fitted circular Gaussian compo- 
nents. During the model-fitting procedure of components, we 
found that inner component B3 (see Fig. 1) was consistently 
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detected in all seven archive data and our own observation (see 
Table 1). At the same time, component B3 had about half as 
much flux density as core when it was ejected from the core. 
Because this was similar to a situation of the innermost com- 
ponent B4, it was continuously detected in at last two obser- 
vations. Component B4 also displays a flux density compara- 
ble to that of the bright core when it is ejected from the core. 
That means we could consider those components B3 and B4 
as real and rule out the possibility that they are just noise on 
the fringes of the strong core component. Table 1 is a summary 
of the model-fitt ing parameters. W e used data from the VSOP 
observations by Lister et al. 200 l l for our component identifi- 
cation. A maximum uncertainty of 15 % in the flux density was 
calculated from the uncertainties of the ampli tude calibration 
and from the formal errors of the model fitting. iFomalont 19891 
presented a method to estimate the position error as Ar - 
where cr is the residual noise of the image after the subtraction 
of the model, is the full width at half maximum (FWHM) 
of the component, and Sp is the peak flux density. However, 
this formula might underestimate the actual position error if 
the peak flux density Sp is very high or the width of com- 
ponent is small. Therefore we derived our position error from 
the formula Ar = Figure 2 presents the core separation 
(the radial distance from the core) as a function of time for the 
individual model components (open symbols denote the data 
at 22 GHz from published literature, solid symbols show our 
model-fitting results i n Table 1 ; stars di splay components from 
the VSOP survey by iLister et al. 200 ll which is used for our 
component identification), where circle, triangle, and square 
symbols denote the Bl, B2, and B3 components, respectively. 

3. Kinematic motion of the jet 

After combining all 14 epochs VLBA observations, we found 
there are two dififerent possibilities to classify the components. 
The first one is related to a simple linear proper motion (see Fig. 
2(b)). The linear proper motion of component B 1 is found to be 
0.21 mas yr consistent with Ijorstad et al. (2001). It is also 
roughly consistent with the result of lKellermann et al. (2004), 
who found a smaller proper motion of 0.15 mas yr~' (4.1 c) 
from a simple linear fit. As for component B3, we found a lin- 
ear proper motion of 0.32 mas yr"' , which is about 50% higher 
than that of component B 1 . Then we found that is always one 
blob ~ 0.5 mas away from the core at diff'erent epochs. If those 
blobs were responding to the same component, it would mean 
that component B2 is quasi-stationary at distances of ~ 0.5 mas 
away from the core, at least during the observation interval. 
However, we note that there is an abrupt transition from the 
year of 1996 September to 1997 June for both components B2 
and B3. Component B2 experiences a sudden rapid fallback to 
the core at an apparent velocity of 0.39 mas yr"' between these 
six months, while B2 is modeled to be in a quasi-stationary mo- 
tion. Meanwhile, component B3 undergoes a fast acceleration 
with 0.78 mas yr"' during the same time interval. This proper 
motion is much higher than that fitted from all observed data 
points of component B3. However, we note that because of the 
rarity of observations, the apparent acceleration when compo- 
nent B3 passes through or gets ahead of component B2, could 



also be a result of the observations' accuracy. More observa- 
tions are needed in the future. 

Although a linear fit can explain the motion of compo- 
nents, there is another possibility, depending on the differ- 
ent classification of the components (see Fig. 2(c)). In this 
case, component B2 seems to remain almost stationary at a 
core separation of ~ 0.5 mas during the period 1994-1997. 
Subsequently component B2 accelerates to 0.20 mas yr"' 
(5.5 c), a velocity comparable to that of component Bl. 
This kind of accelerated motion has sometime been inter- 
preted as a change of the angle between the jet direction 
and the Une of sight ISo ugh et al. 1996t I Vicente et al. 19961: 



Homan etal. 2001 



Wang et al. 200 Ih . Meanwhile, component 



B3 moves with a similar apparent velocity to component Bl be- 
fore reaching a core separation of ~ 0.5 mas. When component 
B3 is ~ 0.5 mas away from the core, its speed is significantly 
low. Thus, it is likely that the two inner jet components (B2 and 
B3) show a similar proper motion mode; both appear station- 
ary at a core-separation of ~ 0.5 mas, but when they leave the 
region of ~ 0.5 mas from the core, they both display velocities 
on the order of 0.2 mas yr ' (5.5 c). Then the linear motion is 
not enough to explain both the core separation evolution and 
the overall path of the compone nts. Therefore we adopt a sim- 
ple helical model presented bv ISteffen et alJ (1995) to fit the 
evolution track of the components. 

In the helical model, the jet can be described by four phys- 
ical parameters, the jet's kinematic energy Eiti„, the opening 
half-angle iff of helix (hereafter called opening angle), the an- 
gular momentum L^, and the momentum along the jet axis P,. 
Keeping three of these four physical parameters constant, four 
possible cases were studied. We found only one case (where 
Ellin, L~, and (/r are conserved) that was physically plausible 
for our kinematic study of B3 1633-1-382. To fully determine 
the curved trajectory with the helical model, 10 free parame- 



ters are required (for an illustration see Fig. 1. in lSteffen et al 
These are: three parameters corresponding to the orientation of 
the coordinates (the viewing angle between the jet axis and the 
observer's line of sight 6; the phase angle of the helix 0o; the 
position angle of the jet axis in the sky plane x), two parame- 
ters describing the initial point of the jet ejection (the ejecting 
time tgj and the distance from the jet axis ro), two parameters 
accounting for the off'set of the observed VLBI-core from the 
ejection of the helix in right ascension (Aor) and declination 
(A(J), and three parameters of conservation (the initial value for 
the Lorentz factor y of a component, the angular velocity wo, 
and the opening angle tfr). The best parameters derived from 
this simple helical fit are listed in Table 2. 

There a re several co nstraints derived by using the same 
methods as ISteffen et al.l (1995): the Lorentz factor (y > 10), 
the opening angle (iff < 5°), the initial radius (ro < 0.1 mas), 
the time of appearance for component B3 (f^y ~ 1995.2), and 
the offset ( ^Aa^ + aS^ < 0.15 mas). We plot the result of 
model fits to the observed trajectory for each component in the 
left panels of Fig. 3. Points with error-bars are observed data, 
whereas the solid lines are from model calculation. 
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The apparent velocity Papp is related to the proper motion 



Papp — 



C{\ +Z) 



(1) 



jPearson & Zensus 1987D . where di is the luminosity distance, 
c is the speed of light, and z is the redshift. Proper motion cal- 
culated from two neighboring data points are assigned to be 
measured at the middle of these two observations. With param- 
eters y and 6 derived from the hehcal model fitting for each jet 
component, we can calculate the predicted apparent velocity 
0anp using 



app 



PsinO 
1 -yScose' 



(2) 



where p - ^JT^-y^ is the velocity of the bulk motion. 0^^^ 
can be compared to the observed Papp, in order to evaluate the 
correctness of this fit. In the right panels of Fig. 3 we show a 
comparison of the apparent velocity P^^p, calculated with the 
helical model, and the apparent velocity Papj, estimated from 
observations. Points with error-bars and solid lines are corre- 
sponding to the observed data and helical model calculation, 
respectively. 

The model parameters of the helical fitting (e.g., y, oj, 
<p) Usted in Table 2 are consistent for the three compo- 
nents. Moreover, the single dish ra dio flux monitoring of 
this source at 22 GHz , reported by (ITeraesranta et al. 1998 : 
Teraesranta et al. 20041) . lasting nearly twenty years from 1981 
to 2000, shows three prominent peaks at about 1986.7 , 
1991.4 and 1995 .3 in its hght curve dBritzen et al. 1999t 
Tiirler et al. 1999h . They can be considered as the ejection time 



tej for each corresponding component, assuming that each new 
component appearance is accompanied by an outburst in the 
total intensity. The ejection times derived by the helical model 
for these three components are 1986.6, 1990.1, and 1994.1. If 
component B 1 moves along a linear path, the t^j derived from 
proper motion extrapolation would be around 1989.7, three 
years later than the flaring time 1986.7. Compared with the 
extrapolated ejection times from a linear fit, it seems that the 
ejection times derived from the helical model are more consis- 
tent. 

After combining the 14 epoch observations of this source, 
we find the kinematical motions of B3 1633+382 can be 
reasonably explained with both a simple linear and heli- 
cal mode. In the linear kinematical motion, the component 
B2 moves at significantly slower apparent velocity than the 
other two components in this compact jet. These low-pattern 
velocity or even stationary features have sometimes been 
explained as standing recoUimation shocks in an initially 
over-pressurized outflow, which have bee n reproduced in nu 



merical simulations of AGN jets (e.g., (IGomez et al. 1995 



Perucho & Marti 20071 iLister et al. 20091) ). Alternatively, it is 
possible that the quasi- stationary motion of the component is 
due to an extremely small viewing angle of the specific ge- 
ometry dAlberdi et al. 20001) . We roughly estimate a Doppler 
factor of approximately several tens by using both the aver- 
age flux density and the average size of the core. Taking the 




Fig. 3. Results of helical model fit to the observed trajec- 
tory and the comparisons of model-calculated apparent veloc- 
ities p^pp for each component (observations shown with black 
points/black solid line with points; model shown with red solid 
Une). 



apparent velocities into account, this implies that the jet is 
aligned at ~ 1° to our line of sight, and the Lorentz factor 
corresponding to the pattern velocity is F ~ 18. Although it 
is an approximate estimation, the direction of the jet from our 
line of sight and the Lorentz factor fitted from linear model 
are consistent with the parameters from the helical model cal- 
culation. As for the helical fitting, on the other hand, it has 
been applied for fitting both the projected trajectory and the 
apparent velocity of components in the compact jet of quasar 
B3 1633+382. The model-calculated projected trajectories also 
well fit the significantly low apparent velocity for components 
B2 and B3 around a core-separation of ~ 0.5 mas. Helical 
motions combined with indications of a time-dependent ejec- 
tion angle can be explained in the context of a binary core. 
Binary black hole (BBH) models have been employed to 
explain similar compon ent motions in other sou rces (e.g., 
1803+784; (Britzen et a L 2010l; iRoland et al. 20081) . 3C 345; 
dLobanov & Roland 20051) ). The component B3 1633+382 
has been shown to exhibit periodicities in the radio 



(IFan et al. 2007n and in the optical (1) lightcurves, while also 
being variable in the y-rays (2). Quasi-periodic variability 
across the spectrum is a further argum ent in favor o f a BBH 
in th e nucleus of this ob j ect (Karou zos et al. 2010t e.g., OJ 
287ilSillanpaa et al. 19881: lUhto & Valtonen 19961 3C 454.3; 
Oian et al. 20071) . Helicity and precession of the jet can be 
explained in terms of Kelvin-He l mholtz instabilities (e.g., 
Camenzing & Krockenberger 19921 Perucho et al. 20061) . pre- 
cession of th e accretion disk (e.g.,I Caproni et al. 2006h . or mag- 
netic torques iLai 20031 Although the stationary proper motion 
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of components in B3 1633+382 can also be explained with both 
linear and helical model in this study, the long-term VLBA 
study has indicated that the standing com ponents may be in 
temporarily quiescent states dLister et al. 2009). Future obser- 
vations will provide new kinematical constraints. 



4. Conclusions 

The quasar B3 1633+382 shows two different kinds of proper 
motion for its components. We note that it could be explained 
as two kinds of kinematical models when we use different clas- 
sifications for the components. In the first case, if component 
B 1 and B3 move outwards with a hnear speed, then component 
B2 stays at a quasi-stationary state ~ 0.5 mas from the core 
for the whole observation interval. Component B 1 has a proper 
motion of 0.21 mas yr"', which is roughly similar to previ- 
ous observations. In the meantime, component B3 also expe- 
rienes a linear expansion with an apparent velocity of 0.32 mas 
yr '. Although this model can explain the motion of the com- 
ponents, there is another possibility. In the second case, com- 
ponent B2 was accelerated after remaining almost stationary 
at a core-separation of about 0.5 mas for several years, while 
component B3 decelerates to a significantly low velocity at the 
same core separation of about 0.5 mas. Then we applied a he- 
lical model under the assumption of three conservative quan- 
tities (the jet's kinematic energy, an gular m omentum, and the 
momentum along the jet axis; Steff en et al.l 1995) to interpret 
the kinematics of the components. The solutions for each com- 
ponent can be used to explain the projected trajectories and the 
apparent velocity with time as well as to demonstrate in partic- 
ular the low apparent velocity of the two components B2 and 
B3. Based on the predictions from the helical model, the com- 
ponents may have the same mode of motion. Future observa- 
tions are needed to provide new kinematical constraints for the 
motion of these three components in this source. 
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Table 1. Gaussian models of VLBI observation. 



Epoch 


Frequency 


Component 


Flux 


R 


6 


Major 




(GHz) 




(Jy) 


(mas) 


(deg) 


(mas) 


1994.17 


22.2 


Core 


0.81 ± 0.12 


0.0 




0.10 






Bl 


0.25 ± 0.04 


0.99 ± 0.05 


-87 


0.48 






B2 


0.30 ± 0.05 


0.46 ± 0.02 


-92 


0.37 


1996.07 


22.2 


Core 


0.74 ±0.11 


0.0 




0.03 






Bl 


0.28 ± 0.04 


1.32 ± 0.05 


-87 


0.78 






B2 


0.08 ± 0.01 


0.53 ± 0.03 


-98 


0.30 






B3 


0.46 ± 0.07 


0.17 ± 0.01 


-86 


0.21 


1997.61 


22.2 


Core 


1.20 ± 0.18 


0.0 




0.07 






Bl 


0.37 ± 0.06 


1.74 ± 0.03 


-87 


0.72 






B2 


0.13 ± 0.02 


0.83 ± 0.06 


-81 


0.67 






B3 


0.13 ± 0.02 


0.39 ± 0.03 


-70 


0.26 


1998.01 


15.4" 


Core 


1.15 ± 0.17 


0.0 




0.11 






Bl 


0.49 ± 0.07 


1.82 ± 0.01 


-85 


0.76 






B2 


0.18 ± 0.03 


0.87 ± 0.02 


-84 


0.67 






B3 


0.14 ± 0.02 


0.50 ± 0.03 


-11 


0.28 


1998.32 


22.2 


Core 


1.38 ± 0.21 


0.0 




0.07 






Bl 


0.35 ± 0.05 


1.90 ± 0.04 


-86 


0.77 






B2 


0.15 ± 0.02 


1.01 ± 0.08 


-78 


0.74 






B3 


0.20 ± 0.03 


0.50 ± 0.04 


-82 


0.33 


1999.04 


22.2 


Core 


1.20 ± 0.18 


0.0 




0.08 






Bl 


0.30 ± 0.04 


2.05 ± 0.04 


-85 


0.91 






B2 


0.20 ± 0.03 


1.18 ± 0.04 


-80 


0.87 






B3 


0.16 ± 0.02 


0.54 ± 0.04 


-82 


0.44 


1999.57 


22.2 


Core 


0.88 ± 0.13 


0.0 




0.04 






Bl 


0.29 ± 0.04 


2.13 ± 0.06 


-84 


0.99 






B2 


0.20 ± 0.03 


1.26 ± 0.03 


-81 


0.63 






B3 


0.14 ± 0.02 


0.42 ± 0.04 


-79 


0.50 






B4 


0.81 ±0.12 


0.13 ± 0.02 


-42 


0.01 


2000.16 


22.2 


Core 


0.82 ± 0.12 


0.0 




0.05 






Bl 


0.25 ± 0.04 


2.15 ± 0.04 


-85 


1.09 






B2 


0.15 ±0.02 


1.33 ± 0.02 


-79 


0.68 






B3 


0.07 ± 0.01 


0.51 ± 0.05 


-79 


0.46 






B4 


0.33 ± 0.05 


0.21 ± 0.02 


-52 


0.12 



Note: a denote the correction for atmospheric opacity for data on 1998 January encountered 
a problem. Thus we used the simultaneously observed 15 GHz data to replace the 22 GHz data. 



Table 2. Parameters of the helical model. 



Component 


tej 


e 


ro 


Wo 


ro 




00 


X 


Aa 


h6 






n 




[^] 


[ly] 


["] 


n 


n 


[mas] 


[mas] 


Bl 


1986.6±0.2 


1.4±0.2 


19.8±1.3 


15.5±0.1 


0.41±0.7 


0.21±0.04 


21±6 


94±2 


0.10 


0.01 


B2 


1990.1±0.2 


1.0±0.1 


19.5±0.4 


23.2±0.1 


0.30±0.1 


0.16±0.06 


65±15 


87.0±2 


-0.05 


-0.04 


B3 


1994.1±0.1 


0.6±0.1 


19.5±0.6 


20.7±0.1 


0.22±0.1 


0.12±0.02 


-79±13 


73±5 


-0.07 


0.02 



Col. (1): observed components; Col. (2): time of ejection; Col. (3): angle between the jet axis and the line of sight; Col. (4): lorentz factor; Col. 
(5): angular velocity; Col. (6): initial radius; Col. (7): opening angle; Col. (8): initial phase angle; Col. (9): projective position angle of jet axis 
in sky plane; Col. (10): offset of ejection in right ascension; Col. (1 1): offset of ejection in declination. 



